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Motivation for This Work 
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NASA 
N+1 Goal 

Open rotors have the 
potential for significant fuel 
burn savings. The challenge 
is to make them acoustically 
competitive too. An essential 
ingredient of such an 
endeavor is the ability to 
robustly predict their noise as 
a function of design changes 
and installation effects. 
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% Fuel Burn Reduction 


Icons represent notional numbers 
based on published information 

2 



Acoustic Prediction Framework 
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Geometry, Flight Path 
& Operating Conditions 



t^u Aerodynamic 
j|p Simulations 




Component Validation Data 
Acquired / To Be Acquired 


NASA Code Toolbox 

ASP Tools FSC APET, RNM RNM, ANOPP SEA/FEM 

TURBO, ADPAC, PAS, 

ASSPIN/PyASSPIN/ASSPIN2, FW-H, 

OVERFLOW, FUN3D 

LINPROP/QPROP, CRPFAN, PAS 


Atmospheric Propagation 
& Terrain Modeling 


Community Noise 
Prediction 


Component 
Noise Source 
Prediction 

/ 

Focus of this presentation 



Airtrame Scattering 
& Shielding Prediction 


Cabin Noise 
Prediction 



Modeling Methodologies 


Noise Sources 

Tone & Broadband 

Thickness (tone only) 

Loading 

Quadrupole 


Note: 

State of the art (or practice) 
for modeling and prediction is 
not the same for all noise 
sources or types. 
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Increasing Complexity 
Increasing Resource Req. 


Computational Aeroacoustics 

Very Pew Attempts to Date 


Acoustic Aerodynamic 

Calculation Step Calculation Step 



A team comprised of NASA GRC and OSU researchers has 
been tackling the aerodynamic analysis of open rotors. 
Currently using TURBO code (NASA & OSU) and Numeca’s 
FINE/Turbo code (NASA) to compute the aerodynamics. 


5 


Aerodynamic Input 




Example: 

TURBO Unsteady RANS Simulation of F31/A31 at Nominal 
Takeoff Condition - Equal RPMs (M tip Hencal = 0.66/0.64) 


12-Bladed Front Rotor 
10-Bladed Rear Rotor 

Measured* 

Predicted 

Thrust 

(Ibf) 

Front Rotor 

303 

304 

Rear Rotor 

305 

309 

Torque 

(ft-lb) 

Front Rotor 

178 

182 

Rear Rotor 

171 

177 

Power 

(hp) 

Front Rotor 

225 

230 

Rear Rotor 

216 

223 


Computed pressure distributions to be compared 
with pressure sensitive paint data. 

Contour plots and integrated quantities shown 
here were computed from a simulation carried 
out at OSU by Trevor Goerig. 


^Preliminary 
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Pressure Loading, Pa 


Computed Blade Loading Time 


Ex.: Front Rotor Time History at an Outboard Point on the Blade 
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Pressure Loading, Pa 


Computed Blade Loading Time Histo 


Ex.: Rear Rotor Time History at an Outboard Point on the Blade 
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Pressure Loading, Pa 


Compute d Blade Loading Time Histo 

Comparison of Time Histories 


95000 



Loading Harmonic Amplitude, dB 


Computed Blade Loading Spectre^ 


There are significant differences between the unsteady loading 
content of the two rotors. Expect to see differences in the 
relative acoustic contributions of the front and rear rows. 


Front Blade Rear Blade 




Loading Harmonic Loading Harmonic 

There is several orders of magnitude difference between the mean and the 
unsteady components. Yet, wind tunnel data indicate that the unsteady 
component can contribute significantly to the overall noise of an open rotor. 


Acoustic Solution - 
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Blade Loading- 


Lighthill Stress 


Blade Normal Velocity 
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Thickness Source 

-T ^(t) (geometric input) Propagation 



F. 


T / \ Loading Source 

-l SjTj (aerodynamic input) 
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Propagation 


T.. 

jk 


/ \ Quadrupole Source 

-l F(tJ (aerodynamic input) 
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Propagation 


Green’s Function 
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Frequency- Domain Solution 

Solution for Isolated Rotor 
(Valid for All Source Types) 

Ref.: E. Envia, AIAA Journal, Vo. 32, No. 2, February 1994 
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Tone 

Frequency 

-i mBQ. t 


Tone Amplitude 


B Blade Count 
Q Rotational Speed 

Airy Function & 
Its Derivative 
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Surface or Volume 
Integral Computed 
Using Quadrature 


Uniform Asymptotic Approximation to Radiation Efficiency Integral (i.e., 
integration over j )• This formula is much more efficient than carrying 
out the integration numerically when the parameter mB is large. 


The representation is valid across the tip speed regime (subsonic to supersonic) and applicable to any 
observer position (near- or far-field). The code based on that solution for thickness & loading sources is 
called LINPROP and that for quadrupole source is called QPROP. The Data-theory comparisons for 
single rotation rotor configurations for both codes can be found in the cited reference. 
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Comparison to Time-Domain 


More Recent Single Rotation Example 
SR7 Propfan Noise Predicted Using LINPROP & ASSPIN Codes 



SR7L: M„ = 0.8, J = 3.26 


SR7L: M„ = 0.8, J = 3.26 
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Loading Noise 

ASSPIN (TK) 

USPROP (TK) 



fOOO 

n o.i ih.i 04 o.y o.ti 0.7 

Fraction of Period 


o v 




Of 0.2 o..< 0 4 f».' f l.o 

Fraction of Period 


o.,s I tv 


/ 


Ext. of LINPROP/QPROP to Open Rot®*s 


OO OO 

* Only loading term , 
needs to be modified p ' M xj 

w=—°° k=—°° 

m is noise harmonic index 
k is loading harmonic index 



Tone Frequency 



Tone Amplitude 


+kB 2 Cl 2 )/ 


Blade counts and 
rotational speeds 
need not be the same 
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mB { - kB 2 
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(mB x -kB 2 ) 

V V ' 

Tone radiation efficiency is mostly 
controlled by this index parameter 
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Weaker loading harmonic amplitude of interaction 
tones relative to that for the primary rotor tones is 
offset by their much higher radiation efficiency. 
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Preliminary Acoustic Predictions 



Measured Sideline Narrowband Acoustic Spectrum at Broadside Position 
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Measured Broadside Spectrum 



2 Measured F31/A31 open rotor acoustic 

spectrum at 60-inch sideline distance. Shown 
levels are not corrected for atmospheric 
| absorption though the effect is small for the 
shaft order range shown here. 
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Tone Extraction from Narrowband Spectrum 
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Extracted tone levels account for the tonal 
so energy spread over 2-3 narrowband bins. 
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Preliminary Acoustic Predictions 


Data-Theory Comparisons Using LINPROP Predictions 


Only thickness and 
loading sources are 
considered in this 
comparison. 
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BPF i~l 


BPF. 


Most strong tones (individual rotor 
blade passing tones and, in 
particular, interaction tones) are 
reasonably well predicted by the 
CFD + LINPROP combination. 
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Preliminary Acoustic Predictions 


Decay rates of measured individual rotor harmonics (i.e., nBPF, & nBBF 2 ) do not conform to the 
expected behavior at subsonic helical tip speeds. Culprit may be either nonlinear propagation 
effects in the nearfield (likely) or contribution from the quadrupole source (not very likely). 
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Shaft Order 


Summary 



An effort has been underway at NASA to assess and improve NASA open 
rotor noise prediction tools. LINPROP and QPROP are among the NASA 
codes for source noise prediction that have been extended and are being 
evaluated for handling open rotor configurations. 


A critical element of the noise prediction process is the computation of the 
unsteady aerodynamic input needed by these codes. 

Preliminary results suggest that LINPROP can predict open rotor interaction 
tone noise reasonably well, but additional improvements may be necessary 
to better match the measured individual rotor harmonics (nBPF! and nBBF 2 ) 
in the nearfield. Given the preponderance of tones in open rotor spectra, 
the highly efficient asymptotic approach incorporated into LINPROP (and 

QPROP) makes it quite suitable for this type of analysis. The bottleneck is the 
CFD input generation process. 


Concurrent efforts at NASA and OSU have been focused on making the 
aerodynamic prediction element more efficient using both fully coupled 
TURBO and “selectively” coupled FINE/Turbo CFD approaches. 

Further analysis and data-theory comparisons are underway to establish the 
accuracy and robustness of the LINPROP and QPROP codes in particular, 
and the FW-H approach in general. 


